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ABSTRACT 



of thesis entitled 

"A Study of the Dynamics Associated with a Non-linear Closed Loop 
System" 

by 

Benedict 0. Olson 

SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE 
OF MASTER OF SCIENCE IN THE DEPARTMENT OF PHYSICS ON MAY 16, 1952 

In the development of computers used in fire control systems, 
there arises the requirement that two vector quantities (represented 
as scaler signals) be added vectorally at right angles. A method that 
was suggested to the author by Dr. Y. J. Liu was to use a closed loop 
system involving the squaring property of electro-magnetic microsyn 
units operating on a shaft. Although the instrumentation of such a 
system is fairly straight forward, the dynamic analysis of the system 
is made complex because of the non-linearities involved. This thesis 
will be in the main a study of the dynamics of the proposed system. 

After a brief discussion of the theory and operation of the 
microsyn units to be incorporated in this system, the performance 
equations along with stability conditions for an ideal system are 
developed in Chapter II. Chapter III then develops a generalized 
approach to the design of a practical system. To keep the dynamic 
qualities of the system constant a square root generator has to be 
incorporated in the feedback loop. What was done in Chapter III was 
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to develop a method specifying how close to an ideal square root 
characteristic we have to be in order to stay within pre-assigned 
tolerances of the natural undamped frequency and the damping ratio 
associated with the system. Chapters IV and V present transient and 
steady state analyses similar to the types used in linear second order 
systems. The final chapter shows the complete design and results of 
a laboratory test system using the general design criteria of the 
previous chapters. 

Thesis supervisor: Dr. Walter Wrigley 

Associate Professor of Aeronautical Engineering 
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CHAPTER I 



Introduction to the Problem 



In the instrumentation of systems related to a wide variety 
of problems, there often arises the need to add vectorally two quan- 
tities at right angles. In general, the two vector quantities are 
represented as scalar signals, and it is the requirement of some 
instrument or system to treat these signals as vectors and perform 
the right angle addition. In this thesis, we shall represent the 
two signals as voltage quantities £^nd Ey. The desired result will 

r 2 2 

then be Eq = J Ex ^ where again Eq is a scalar quantity. It is 
to be noted here that we are not interested in the actual position with 
respect to inertial space of the two vectors, but merely in restoring 
their right angle relationship with respect to each other. 

The quantities E^^ and Ey will fluctuate either periodically 
or haphazardly, but they will never change direction. The reason why 
and Ey can only be fluctuating direct voltages of the same sign will 
be shown later in the chapter. A major portion of this thesis will be 
devoted to a study and analysis of the dynamics associated with the 
proposed system. In short, we shall attempt to answer the two questions, 
*^how does the system respond to fluctuating inputs?"- and given certain 
tolerances on how well the system should perform, "how can we design a 
practical system?" 

One method in use at present to obtain the result 

I 2 ^~2 

Eq “ ^ E^ ^ ^ is to use a resolver nulling servo combination. The 
system is shown schematically in Fig. 1-2, 
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DESIRED RESULT : 
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FIGURE l-l 

RIGHT ANGLE RELATIONSHIPS OF QUANTITIES TO BE 
COMPUTED 




USUAL METHOD OF SOLVING PROBLEM 
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and Ey are root mean square values of alternating sinus- 
oidal voltages appearing across the stator windings I and II. 

Eq 2 voltages appearing across the rotor windings as shown. The 
axes of the two stator windings are located 90 mechanical degrees apart 
as are the axes of the two rotor windings. Because of the sinusoidal 
relationship of flux density in the air gap as a function of angular 
position of the stator periphery, the following relations hold for 
open circuit rotor windings. 

Eo 2 " Ex a *4. Ey sin a 1-2 

In the case of the resolver nulling servo combination as 
shown in Fig. 1-fe, the shaft of the resolver is driven in a direction 
by the servo motor so as to null EqI* nulling servo amplifier 

has a large enough gain, EqJ can, to a first approximation, be con- 
sidered zero but it is understood that Eq^ = 0 only if the gain of 
the nulling servo amplifier is infinite. Thus equations 1-1 and 1-2 
become 

Eo2 “ Ex cos a + Ey sin a 1-2 

0 = - Ex OC + Ey cos Ct 1-3 

Squaring equation 1-2 and 1-3 and adding, we obtain 

and thus we have our desired result 

^ - h -- 

This system as snown in Fig. 1-2 is a true servomechanism 
and its characteristics will be used later as a comparison showing the 
differences between this true servomechanism and the closed loop system 
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to be described in this thesis. 

We shall hereby name the proposed system the orthogonal vec- 
tor component sunBiing system. A simplified diagram is shown in Fig. 1-3# 

At this point, we shall neglect the damping and inertia terms 
and merely show how the steady state iresult is obtained. In Fig. 1-3 
the two inputs and ly are acted upon by an operating component (oc) 
such that an output torque is produced on a shaft (torque summing 
member) which is proportional to the square of the input signals. The 
shaft rotates through an angle to which is attached an operating 

component such that an output voltage is produced proportional to the 
angle of rotation* The square root of this voltage is taken and fed 
to the shaft through a third torque producing device which depends 
on the square of its actuating input quantity* 

In the steady state, that is, assuming all transients have 
died out, we have the three torques and acting on the shaft* 

-h a Mj 1-5 

2 * 

S(oc)l [e2,M] + Hoc)2 [e2,n] =* Hoo)3 [ehvi] 

Assuming the sensitivities giving torque out for voltage 
square in are equal, we arrive at the desired result, namely, 

Eq2 « E x2 + e/ 

^ + V 

The square root function in the feedback loop cascades with 
the squaring device to produce a linear device in the feedback 

* 

Refer to the i^pendix for explanation of terminology used. 
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FIGURE 1-3 

PROPOSED SYSTEM IN SIMPLIFIED FORM 
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loop. Equation 1-6 shows the very important result that in the steady 
state, the correct answer to the problem will be given regardless of 
the nature of the device used in the feedback loop. However, in the 
next chapter where the complete performance equations of the system 
will be developed, we shall show that a perfect square root function 
is required in the feedback loop if the dynamic qualities of the system 
are to remain constant. The remainder of this chapter will be devoted 
to a description and theory of the torque producing generator and the 
angle input voltage output component commonly called the signal 
generator. 

Both the torque generator an^ signal generator belong to a 

1 ^ 

family of devices knovm as microsyn electromagnetic units. These 
units consist of a laminated ferromagnetic rotor of 180^ symmetry with 
two poles and a laminated ferromagnetic stator of 90^ symmetry having 
four poles about which are wound various arrangements of electrical 
coils depending upon the application of the unit. For the types of 
microsyns used in connection with this thesis, the rotor is allowed to 
rotate - 10^ from a zero rotor position, which is geometrically defined 
when a rotor pole is half way between two adjacent stator poles. Fig. 
2-4 shows the general construction of a microsyn. torcue generator. 

The torque generator produces a torque w^hich, in the ideal 
case, is proportional to the square of the current through one axis 
and is independent of rotor angular position over the useful range. 

«(tg) = 1-a 

V/e shall attempt to clarify this relationship from energy 
considerations in a magnetic field. 

As a starting point, we note the fact that the total 
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SHAFT ON WHICH ROTOR IS MOUNTED 
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electrical energy into the system equals the magnetically stored 
energy plus the work done by the system, neglecting heating losses. 
That is, 



'“■'(.lec) = 

The instantaneous electric power delivered to the torque 
producing coil is ei. Thus, after a time t, the electric energy 
into the system, is 



According to 



(elec) 

Faraday’s law 
^(ind) = 




ei dt 



1-10 



where 



Hence 



Uliere 



®(ind) 

N 



voltage induced by motion 

the number of conductors 

magnetic flux cut by the conductors 



(elec) 




1-11 



is flux at t = 0 
02 is flux at t = t^^ 

We shall go from 01 to 02 paths as shown in 

Fig. 1-5. First, following path 1, the current is maintained constant 
and hence 



(Eft) (elec) = 

This will correspond to a certain rotation of the rotor 
to a point corresponding to the flux 02* following path 2, 
we let the current drop to zero, and if we neglect any hyster6ses 
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d) ► 



FIGURE 1-5 

AMPERES VS. FLUX IN A MAGNETIC MEDIUM 
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effects, the flux too will drop to zero. At this point the rotor is 
rotated through the same angle as in the first case. However, no work 
is done because there is no force. Then the current is increased until 
the point (02> reached. In the second case, the magnetic energy 

can be written in the following form. 




1-13 



Going from 0 to zero, denote ^ where is the 



slope of this line. Similarly, going from 0 to 02? <^enote ^2 = ^ 
the slope of this straight path. Then 1-13 becomes 




1-14 



But 




Therefore 






or finally 



(ai)(eleo) = 



1-15 
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Equation 1-12 represents the area under path 1 while 1-15 



represents the difference between the areas under the two sections of 
path 2. 



that the mechanical work done and the magnetic energy change are equal 
in magnitude for a given amount of input electric energy. 



field equals the work done under the conditions previously described 
enables us to calculate the work done by finding its equivalent in 
terms of the magnetic energy storage. Then, the torque can be found 
by taking the derivative of this amount of energy with respect to the 
rotor angle. 



assumptions : 

1. The reluctance of the magnetic circuit is concentrated almost 
entirely (over 90^) in the air gaps. As a consequence of this assumption, 
we shall assume that the entire magnetic drop around the magnetic loop 
produced by the sources of magnetomotive forces at opposite poles occur 
essentially across the air gap. Thus, we can assume to a good approxi- 
mation that all the magnetic energy is stored in the air gap. 

2. Leakage flux will be neglected. 

From reference (2), the magnetic energy stored in a given 
volume is given by 



where V is the volume, B is the flux density. As mentioned above, 
this energy will be assumed to be stored in the air gap. 

Hence: 



Since 1-15 is exactly half that found in 1-12, we conclude 



The fact that the magnetic energy stored in the magnetic 



Before proceeding further, we shall make the following 




HdB 



1-16 
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V 



1-17 






— f 

® B. 



B, 



BdB = 



8 ^iUo 



(b/ - 



where 



is the permitivity of air in MKS units. 



(En) 



LA 



(mag) 87T ii 






1-18 



where L is the Histance between the rotor and stator. 

A is total area between the rotor and the two opposite 
stators. Since 



0 



(mag) 



mmf ^ nimf 

reluctance ^(raag) 

^ (change in nnnf)^ 

8 7T A R (mag) 

^qA U(mag)^ 

8 7TL 



where U(mag) is the magnetic drop across the air gaps. This drop 
equals to a good approximation the impressed magnetomotive forces. 
Since mmf = 4^^Ni, we have finally 



^^)(mag) 



2 ^ M aN^ i^ 
L 



W lT=no. turns 1-20 



The torque M then is ^ where A(i.ot) is ^he rotor angle. 

d A(rot) 



Thus it is shown that the torque is proportional to the 
current squared. An important consequence of 1-20 is the torque is 
independent of the direction of the current. Thus altei*nating current 
can be used as well as direct current to provide torques. This also 
shows why only absolute values of and Ey can be considered in the 
orthogonal vector summing system. 
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The change of area A has a linear relationship to a change 



in rotor angle 1-20 becomes 

, , - 2JTflo»^A A(rot> 

L 



where K can be calculated from the dimensions of the particular 
microsyn. 

Finally, 



M, , - M = 

^'(tg) - aA(rot) L 



^(tg) 



.2 

1 



1-21 



where we shall define 

27r^oKN^ 

^(tg)[i*;M] ■ L 

The signal generator provides an open circuit voltage which 
is proportional to the amount of angular rotation of the rotor from 
its null position. Thus 

®(sg) = 2(sg)[A;e] 

Fig. 1-6 shows the details of the signal generator. There 
are two windings per stator pole an*^ these windings are connected in 
two series circuits, each circuit using a similar set of windings. 

The primary circuit is excited by an alternating current which sets 
up flux patterns as shown. If the previously mentioned torque generator 
has equal excitation currents on both axes, the flux pattern will be 
the same as that set up by the primary excitation current of the sig- 
nal generator and thus the resultant torque on the signal generator 
is zero since the torques on the tv;o axes are equal and opposite. There 
will be no torque contribution from the secondary winding since this 
vrinding will carry no current. The series windings of the secondary 
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A C 

INPUT 

SIGNAL 




NOTE-- 

1. HEAVY ARROW ON EACH POLE INDICATES DIRECTION OF 
POLE FLUX WHEN INSTANTAN EOUS DIRECTION OF 

IS THAT SHOWN. 

2. ALL INPUT COIL WINDINGS HAVE THE SAME NUMBER OF 
TURNS AND ALL OUTPUT COIL WINDINGS HAVE THE SAME 
NUMBER OF TURNS- 



FIGURE 1-6 

HICROSYN SIGNAL GENERATOR 
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coils are so arranged that the instantaneous induced voltages in these 
secondaries are added in opposite phase for adjacent poles. 

Thus, 

®(out)(gg) = + %) - (^2 1-23 

where is the rms induced voltage in coil n. n = 1, 2, 3, or 

\\laen the rotor is in its null position (i.e. midv;ay between 
two adjacent stator poles), we see from 1-23 that the signal generator 
output voltage is zero because all the induced secondary voltages are 
zero due to the symmetry of the system. 

A positive rota.tion of the rotor will be defined as a 
rotation tending to increase the overlap area betv;een the rotor and 
stator poles 2 and 4 and correspondingly, to decrease the overall over- 
lap between the rotor and poles 1 and 3, as shown on Fig. Thus, 

for a positive rotation, the reluctance of the magnetic circuit inclu- 
ding poles 1 and 3 increases while that of 2 and 4 decreases. Hence, 
since the magnetomotive force at all the poles is constant, the flux 
in poles 1 and 3 decreases while that of 2 and 4 increases. It is to 
be noted that the change in reluctance and therefore the change in flux 
is a linear function of rotor angle. The induced root mean square 
voltages are proportional to the flux for constant excitation frequency. 
Thus for a positive rotation, the induced voltages in poles 2 and 4 
are the larger and the output voltage has a phase similar to the phase 
of the voltages in poles 2 and 4* For negative rotations, is 

of opposite phase and of magnitude given by 1-23. 

The flux in the poles of the signal generator have very 
nearly the same phase as that of the alternating excitction current, 
the only difference in phase being that attributable to the core losses. 
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This flux is directly proportional to the number of primary coil turns, 
and the magnitude of the excitation current and inversely proportional 
to the circuit reluctance. The secondary induced voltage has a 90^ 
phase shift with respect to the flux and its value is proportional to 
the magnitude of the flux, number of secondary turns, and excitation 
frequency. We can thus define the output voltage as 

®(out) “ ^(sg) [A;i,n;e] ^(ex) ^(ex) ^(rot) ^ ^ 

= S(sg)|A;e] A(rot) 

where 

S(sg) [A;i,n;e] ^(ex) '^(ex) = ^(sg) [A;eJ 

when the excitation current and frequency are fixed, and ^{rot) is the 
angle of the rotor measured from the null position with a positive 
rotation having been already defined as a rotation which increases the 
area under poles 2 and 4. 
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CHAPTER II 



Proposed Ideal System 



The orthogonal vector summing system to be described in this 
chapter will be considered ideal for the following reasons: First, the 

square root generator will be considered as producing an exact square 
root over the entire operating range. Second, all uncertainty torques 
such as those attributable to friction, hysteresis, leakage flux, plus 
interfering torques such as the acceleration of the case with respect 
to inertial space will be neglected. 

Fig. 2-1 shows a line schematic diagram of the proposed 
assembly and Fig. 2-2 shows a functional diagram of the orthogonal 
vector summing system. Since this thesis is intended to be primarily 
a study of the dynamics associated with the proposed system, the des- 
cription of the actual assembly will be mainly to establish a clear 
picture of what physical processes are actually taking place. The 
analysis to follow will be largely mathematical, but nevertheless the 
physical interpretation of the derived performance equations will be 
continually stressed. 

The assembly shown in line schematic form in Fig. 2-1 is 
mounted inside a sealed cyclindrical container and designed to operate 
in an upright position. The torque summing shaft contains the rotors 
of the microsyn units as well as the rotor of the viscous damper. The 
fluid in addition to providing a damping torque, also serves to reduce 
friction torques at the pivots by buoying up the shaft by a sufficient 
amount to remove the weight of the shaft from the pivots. A float is 
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LINE SCHEMATIC DIAGRAM OF THE ORTHOGONAL VECTOR SUMMING 

SYSTEM 



{ tsm) 
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FIGURE 2-2 

FUNCTIONAL DIAGRAM OF IDEAL ORTHOGONAL VECTOR SUMMING SYSTEM 



clamped to the shaft to provide the necessary amourrt of lift. The 
damper consists of two concentric cylindrical surfaces one of which 
rotates with the shaft while the other surface is fixed to the case. 

A damper adjustment is provided which changes the amount of surface 
between the concentric cylinders. Because the viscosity of the damping 
fluid is critically dependent upon temperature, the unit is maintained 
at a constant temperature plus or minus a certain tolerance by means of 
electrically heated coils. A thermostat is provided which sets up a 
cycling process in the application of electric power to the heating 
coils. A portion of Chapter V will be devoted to an analysis of how 
changes in damping effects the dynamics of the system, and thus 
determine criteria for specifying how well the temperature must be 
regulated. 

We shall now develop in detail the performance equations for 
the ideal system referring to Fig. 2-2. Because of the fact that we 
shall regard the input signals as voltages, a current source has been 
provided to drive each of the torque generators. 

By Newton^ s second law, the product of the moment of inertia 
of the shaft plus rotors times the angular acceleration of the shaft 
equals the sum of the torques applied to the shaft. 

ll-(ca)] (OA) * *(rot)) = ^^(vis) A(rot) ♦ ^(tg)2 

♦ "(tg)3 2-1 

Where I = moment of inertia of shaft including rotors of the daji^er 
and microsyn units. 

A(rot) * angle through which a rotor has turned from its null 
position. We shall assume all rotors are lined up to 
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have the same null angle with respect to case coordinates, 
A positive rotation of has been defined as that 

direction tending to increase the overlap area between 
the rotor and stator poles 2 and 4. 

^[l-(ca)j (OA) “ component of angular acceleration of the case with 

respect to inertial space about the output axis, that 
is^ about the direction indicated by a line parallel 
to the shaft. 

= The damping coefficient of the system consisting of 
the following two parts: Cd(fixed)> damping co- 

efficient of the fixed section of the viscous damper 
and C(5(var) variable section of the viscous 
damper. 

Torque generator 1 and torque generator 2 will be actuated 
by the input signals and will be defined as producing positive torques, 
i.e. torques tending to produce a positive rotation of shaft angle. 

On the other hand, torque generator 3 will be actuated by the output 
of the system and will be designed to produce a negative torque. 

Rearranging 2-1, 

I A(rot) ^ Hrot) " ^(tg )3 = ^(tg)i ^(tg)2 " ^ '^[l-(ca)] (OA) ^ ^ 

From chapter I, we have defined 

M(tg) = S(tg)[i2;M] 2-3 

Define S(cs)|^e;i] e’ current source 

for voltage input current output. Hence 

^(tg) = ^(tg)[i^;M] ^^(cs)[e;i] ^ 2-4 
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E. 



re root 



Define , sensitivity of squs 

1 

generator for (voltage)^ in, voltage out. 

From Chapter I, we have defined the signal generator sensi- 
tivityas S(3g)|-A.ej= -f where 



’(sg)[A;e] ^(sg) [A;i,n,e] ^(ex) "(ex) 



where i(gx) excitation current of the signal generator primary 

winding and is the frequency of the excitation current. 

Using the above definitions for the parameters of the system, 
2-2 becomes 



I A(rot) ^ Crf ^{rot) 





S(tg)i[i^;M] S^(cs)i[3;i] ^x^ 2(tg)2fi^;Mf (cs)2[e;i | ^ ^ 



- I A[l-(ca)] 

Since 



2-5 



^(tg)i “ ^(tg)2 ■ ^(tg)3 = 2(tg) 

and 

S(cs)3 = ^(cs)2 ■ ^(05)3 ■ ^(cs)’ becomes 

I A(rot) + ^d ^(rot) ^(tg)j^i^;Mj ^ (cs)|^e;ij ^o 

^(tg) [i^;M]^ (cs)[ejij(^ ^ ) 

Also, 

£0^ : ^^(sr) je^jej ^(sg) [A;e] ^(rot) 
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Using 2-7 and defining 

2(ovss) [A;M] “ ^(tg) |^i^,M] ^ (cs)[e;ij ^ (sr)[e^;e] ^(sg)[A;e] 

the sensitivity of the orthogonal vector summing system for shaft 
angle input torque output, 

2-6 becomes finally 

2 

^ *^*(rot) ^ ^(rot) ^(ovss) |^A;Mj (rot) (tg)[^i ;mJ (cs)[^e;ij 

I Hl-(ca)] <“*> 

2-9 



In the discussion of the ideal system to follow, we shall 
assume that the acceleration of the case with respect to inertial 
space is zero. This acceleration produces an uncertainty torque 
which will alter the desired result. Thus, equation 2-9 simplifies to 
I A(rot) + Hrot) * ^(ovss) [A;m] ^(rot) = ^(tg) ^ (cs)[e;ij 

* V) 

2 

Relating and and the derivatives, we obtain 

Eo^ 

" S^(sr)[e2;e] ^(sg) [ft-e] 

. _ Eq^ 

^ (sr)[ez;e] ^(sg) [Aje] 

• • 

Eo^ 

where the independence of the sensitivities of the signal generator 
and square root generator with time are assume for the range of fre- 
quencies to be considered in this thesis. This statement will be 
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clarified later. 



Using the above expressions for ^(^^ot) derivatives, 



I Eq^ + ^ Eq^ + ^{ovss) [ajM] 

dt^ '■ ** 



2 _ 



dt 



^(ovss)[a|m] 



or 



En 



Cd 



^ E^ 



S(ovss) pl;Mj dt^ 



^(ovss) (A;Mj 



IT ^ V ‘ 

= Ex Ey 



2-12 



2 2-13 



Equation 2-13 is a linear differential equation of the 
second order with constant coefficients in the ideal case. Using 
the standard form of the second order linear equation, we have 

1 d^ £ 2 ^ 2 (DR)(ovs3) 1. Eo^ + 4 Ey^ 2-14 



dt 



Wr, 



dt 



where Wj^ = the natural radian frequency of the system. 

(DR)^Q^gg^ - the ratio of the amount of damping present in 

the system to that amount required for critical 
damping. 

^DH^and Wy^^when knov^n for a second order system, completely determine 
its djnamic behavior. Comparing 2-13 and 2-14, we obtain, 

^(ovss) [aim] 



or 



J S(ovss) tf 



2-15 



or 



2(DR) (ov3s) _ ^d 

^(ovss) ^A|M] 



Cd I -L 

(dr)(ov3s) = t V TsTZ : 



^(ovss) [A|Mj 



2-16 
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Since the two input quantities and Ey react upon the 

shaft through identical operating components, we can let 
P 2 2 

V ^ - Ein in further discussion of the dynamics of the system 
without any loss in generality. Thus, 2-13 becomes 



I ^ E ^ Cri d_ 

^(ovss)p;Mj ^(ovss) 

or its equivalent form 





1 


d2 






Where 


Efn i 
I 


^(ovss 


) [AjM 


or 





Eo" 



= E, 



Wn 



dt 



In 



2-17 



2-18 



2-19 



="(=r)[ei;e] "(=s) [A;e] ‘(rot) = 'tc" * 



2-20 



From the steady state value of 2-13 we obtain 

Eo^ = E^2 ^ Ey2 = 2-21 

2 

Thus, it can be seen that Eq and hence Eq, is not 
critically dependent upon the value of If the product, 

S^(sr)[e5;e] ^(sg) [A;e] cJ’-anged, A(j.ot) merely readjusts itself, 
within its constraints of 10^, to a new steady state value corres- 
ponding to the required steady state, value of Eq • However, it is 
very important to note at this point that a change in the product, 
^^(sr)[ei;ej ^(sg) [A;e] > alter the steady state conditions, 

the dynamics of the system are definitely altered as can be seen from 
2-13, 2-15 and 2-16. The fact that a perfect square root function can- 
not be generated means that the sensitivity product will be continually 
varying over the operating range and the next chapter will go into the 
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details of the relationship between the dynamics of the system and 
the nonlinear device. It is to be noted here that the torque generator 
can be regarded as an aLnost perfect squaring device for the range of 
actuating currents used here. That is, the B-H curve is linear for 
the range of currents used, (O to 40 ma) but there are small uncer- 
tainty torques due mainly to hysterises and to reaction torques. 
Reaction torques are present when the torque produced is not completely 
independent of rotor angle. For the present, uncertainty torques will 
be neglected. 

When the sensitivity product is 

made very high, decreases to a point where it becomes essentially 

an error signal analogous to that of a conventional servomechanism, 
similar to that shown in Chapter I. The system will be analyzed and 
designed, however, from the viewpoint that made to cover 

the full 10^. This condition is what makes this system unique from 
the conventional servo systems. 

A very important point on the stability of the system must 
now be brought up. At the null position, a condition of instability 
exists such that if goes the least bit negative, the system 

drives itself into the stops provided just above the 10® mark of 
rotation. 

Fig. 2-4 shows a stable condition where the rotor has 
rotated through a positive angle. The rotor has rotated to a position 
so that the signal which it produces causes a current iQ to be equi- 
valent to iin ^nd thus the positive torque produced by iin equals the 
negative torque produced by i^. If by some disturbance goes to 

^(rot) ^ ^(rot) ^in same, then io goes to io + d io and 

produces a torque tending to pull + d back to its initial 
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FIGURE 2-4 

TORQUE GENERATOR IN STABLE REGION 
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value of ^{rot)* ^(rot) disturbed to a position - d ^(rot) 

lln constant, then io becomes i^ - d i^ and the torque produced in the 
1, 3 axis is reduced thus allowing - d A(j.Qt) rotate back to 

A(rot)* This again clearly shows the stable condition of Fig. 2-4. 

The configuration shov^Ti by Fig. 2-4 corresponds to a potential minimum 
which, of course, will thus yield a stable system. 

Fig. 2-5 shows the condition of instability. Assume, initially 
the conditions of Fig. 2-5. produces a signal which will pro- 
duce the same current, io, s-s + ^{rot) * condition as shown can 

exist just like a pendulum at the top of its arc. However, if 

- A(j,Qt) disturbed to a position, ^ ^(rot)^ ” 

- A(rot) - ^ A(rot)^ goes to io + d io thus increasing the negative 

torque which in turn increases - ^{rot) system will rotate 

negatively until finally hitting the stops. On the other hand, if 

- A(rot) is disturbed to - (A(rot) " d ~ ^(rot) ^ ^(rot)» 

io goes to io - d i© and the negative torque produced by Iq is 
decreased thus tending to further reduce - ^{rot) null is 

reached. A positive angle is built up on the stable side until finally 
the condition of Fig. 2-4 is reached. Here, we have represented a 
potential maximum condition, it being similar to a physical pendulum 

at the top of its arc. The slightest disturbance in either direction 
results in drastic changes in the motion of both the given computing 
system and the pendulum. 

As a third case, assume initially there is no input excitation 
and = 0. Next assume builds up linearly from zero. As 

builds up from 0 to d iin> there is a positive torque produced by the 
current flov;ing in coils 2 and 4 (see Fig. 2-4). At the same time, 
there is no current io in coils 1 and 3 ^nd hence no negative torque 
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FIGURE 2-5 

THE UNSTABLE REGION OF SHAFT ROTATION 
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is produced. This assures a positive rotation and a stable condition. 

Finally, consider the case where = 0 and = 0. 

If A(rot) disturbed to a position - d there results a signal 

producing an output current d Iq ^.nd a negative torque which will in- 
crease the negative angle. Thus, we have the unstable case where ^(rot) 
progressively increases until the stops are reached. 

We conclude from the above analysis that a constraint has 
to be applied to the shaft so that negative angles as shown in Fig. 2-4 
cannot be built up. If the constraint in the form of a mechanical 
stop is placed exactly at the null, then the system will produce an 
Eq for Ex and Ey even in the limit as E^ and Ey both approach zero. 

If the stop, however, is placed at a point allowing a small negative 
angle - d then the chances are that initially the rotor will 

have passed into this negative region in the absence of input signals 
as was pointed out in the last case of instability. Thus, it will 
take a positive torque produced by the input signals just barely 
greater than the negative torque produced by - d in the feed- 

back loop to rotate the shaft into the stable region. It is difficult 
physically to put a stop exactly at the null, and in the lasA chapter 
it will be shown why a stop can be put on the positive side of the 
shaft rotation and still be able to solve the orthogonal vector addition 
even as the input signals approach zero. This method will be feasible, 
however, only if the swing of ^(rot) fairly large in comparison to 
the angular position of the constraint, and also a certain amount of 
tolerance is allowed on the dynamical characteristics of the system. 

Equation 2-18, while linear when expressed in terms of the 
squares of the input and output signals, is nonlinear when expressed 
in terms of the input and output signals themselves. 
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An insight into the complexity of the nonlinear form of 



2-18 can be gained if it is rewritten as follows. 



dt 




2Eo 



dEp 

dt 



dt2 



Eo^ = 



2E. 



d^Eo 

dt2 



2(^)^ 

dt 



Hence 2-18 becomes, upon rearranging, 

2 2 p 

Eq ^^Eq ^ (^)2 ^ 2(DR) Wn £3^2 + E©^ = E^^ (t) 2-22 

dt2 ''dt ' dt 2 2 

It is to be emphasized here that Wj^ and (DR) refer only to 

the linear second order system defining the square of the actual input 

p , 

and output quantities. Since the linear equation defining Ep is 
convenient to use, practically all the information to be obtained 
from the system will be in terms of the squared functions. The last 
step will then be to analyze the square root of the response function 
associated with the linear second order equation. 
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CHAPTER III 



The Square Root Generator Problem 



The performance equations of Chapter II were developed 
assuming an ideal square root generator in the feedback loop. As we 
have already mentioned, it is very difficult to find in nature a phy- 
sical device that will generate such a curve as shown in Fig. 3-1, 
although of course, there do exist many non-linear de.vices or means 
of synthesizing such devices. 

Equation 2-8 defines the overall sensitivity of the system 



as 

^(ovss) [A;M] “ ^(tg) [i^;M] ^ (cs)[e;i] ^ (sr)[e 2 ;ej (sg)|^A;ej 

Since the torque generator follows a rather precise square 
curve for the range of currents used, only a perfect square root function 
will maintain j constant over the entire range. Physically, 

^(ovss)[a-m] analogous to the spring constant of a mechanical 

harmonic oscillator. This is shown clearly in the equations for the 
natural frequency and damping ratio where 



(ovss) 



S(ovss) [A;m] 

f 



ay ^ 

2 ^ ^(ovss) [ajMJ 



where and (DR)(Q^gg) refer to the performance equation relating the 
squares of the input and output functions. It is to be noted again here 
that the static output of the system is not dependent upon j. 

If instead of a square root function we had a linear function, S(ovss)[A;M 
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(VOLTS) 



note: 

NUMERICAL CONSTANTS FOR 
THIS FUNCTION HAVE BEEN 
DERIVED IN CHAPTER IZI. 



70.7 V 




FIGURE 3-1 

IDEAL NON-LINEAR CHARACTERISTICS IN FEEDBACK LOOP 
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would correspond to a very stiff spring for large rotor angles, and 
hence the Wj^ of the system would be high, and the (^R)(ovss) 
low. For low rotor angles, the reverse conditions would be true* 

V/e thus have to resign ourselves to the fact that in an actual 
physical system Wj^ and (DR)(^^gg^ are not going to remain constant* 
and (l^R)(ovss) corresponding to the ideal system will now be defined 
as a reference natural frequency and a reference damping ratio, 
abbreviated, 

^n(ref) (ovss) (ref ) 

In attempting to solve the square root generator problem, 
we shall apply a method of perturbations to a simplified system and then 
critically examine the conclusions obtained from such a system. 

The simplified form of the perturbation theory will be a means 
of linearizing the non-linear equation relating Eq ^nd ^±r\* The pro- 
cedure will be as follows. We shall assume an arbitrary static operating 
point. Hie input signal will then be varied by a small amount and the 
corresponding small change in the output determined. Expansions by 
Taylor’s series will be used when necessary and because the variations 
will be very small, ^ *s of higher order than one will be neglected. 
i/vTien the so-called perturbed equation is subtracted from the original 
equation, a linear equation results. 

Rig* 3*2 shows the orthogonal vector summing system in some- 
what simplified form where instead of square root and signal generators 
we have merely specified a device w-hose output, Eq is a function of the 
rotor angle, A. 

From Fig. 3-2, and using Kevvton’s second law of motion 

I A = - Cd A ^ Ml * M2 3-1 
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(cs),i:e;U !ill! J S(tg),n?Ml 



( 
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THE NOH-IDEAL SYSTEM 



or 



I A * Crt A - M2 = Mx 3-2 

M2 = - 2(tg)[i2;M] 3^(cs) e;i] 

Ml = ^(tg)[i^;M] ^^(cs)[e;i] ^in 



We shall abbreviate the shaft angle A(rot) ^ this 
chapter because of the large number of times that it will be used. 

Using the fact that 

Eq = f(A), 3-2 becomes 

I A + Cd A . S(tg) s2(^3) [f(A)] 2 3 s(tg) S^(cs) ^^in) 3-3 

where we have assumed that the torque generator sensitivities and 
current source sensitivities are equal. Using the method of pertur- 
bations, we let assumne a new value ^ ^in ^here A E^n ^ 

small change in the static value of the input signal. Then, corres- 
ponding to this change the output rotor changes to A + A A and the new 
output can be written f(A -i- A A) . 

Equation 3-3 for this change in input then becomes 



I(A . A A) ^ C(j(A + A a; 

Expanding f(A 
series we obtain, 

f(A . 2\A) = 



^ S(tg) S^(cs) = ^(tg) ^^(cs)Pin-^^^in 

3-4 

■ ^ a) in terms of f(A) by means of Taylor's 



f(A) . 



d f(A) 
d A 



A A * 



3-5 



Thus, 



where 



We have neglected 4 's of higher order than one. 



[f(A + Ak)]^ = [f(A)]^ * 2f(A) f'(A)AA + 



f(A) = 

d A 



3-6 
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and 



(Ei„ * 4 . 2 4 E 



in 



3-7 



Substituting 3-6 and 3-7 into 3-^ gives 

2 



I A + i4a + Cj A + Cjj/lA + S(^g) ^^(cs) ■'■^(tg) ^ (cs) 



S(tg) S^(cs) E^in 2 S(^g) 3-8 



Subtracting the perturbed equation, 3-^ from 3*“3, the 
following linear equation with constant coefficients is obtained for 
every Ej^j^ and its corresponding rotor angle A. 



I 4a > Cd^A ^ 2 S(tg) s2(^g) f(A) f'(A)4A = 2 S(tg) s2(^g) E.^ 4 E-^^ 

3-9 



Rearranging 3-9 into its conventional form and comparing 
the coefficients with those of the standard second order liner differ- 
ential equation as was done in Chapter II, give the following two 
equations. 



2S 



I 

(tg) ^ (cs) 



f(A) f'(A) 



4A 



2S 



Cd 

(tg) ^^(cs) 



TTa) f'(A) 



.. 4a +4a 



2in 4Ein 
“ f(A) f'(A) 



3-10 



1 4 A 4 2(DR) 4A4 4a = ^ %n 

^ Wn f(A) f(A) 



2 S(tg) S^cs) ^(A) f'(A ) 
I 




and 



3-11 



3-12 
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(^^)(0V3S) 



£d 

2 



/ ~ 

2 S(tg) S2(cs) f(A) f'(A) I 



3-13 



The fact that f(A) and f*(A) occur as products will be very 
important in determining tolerances on the dynamic constants of the 
system. 

We shall now check the above analysis for the ideal system 
described in Chapter II. 

V/e have shown in Chapter II that 



■ S^(sr)[e2;a] ^(sg)[A;e] ^ 

Hence 

Eo - S(sr)[e5;e] ^^(sg) [A;e] 

Thus, in the irleal case, we see that 

f(A) - ^^(sg) [A;e] 

and 

f'(A) = i S(si.)[^el;e] ^^(sg)[A;e] ^ ^ 

and finally, we have 

f(A) f‘(A) = k ^^(sr) [ei;©] ^(sg)[A;e] 



3-14 

3-15 

3-16 

3-17 

3-18 



which, when substituted into 3-12 and 3-13 gives 

S^(cs) S^(sr) ^(sg) ~ 3 J ^(ovss)[a;mJ 



3-19 



and 



, C(} 

(DR)(ovss) - — 



I (S(tg) ^^(cs) ^^(sr) ^(sg)^ 






(S(ovss) ^A jKJ^ 



3-20 
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Equation 3-19 and 3-20 are identical with those developed 
in Chapter II for the ideal system using the conventional approach. 

As mentioned before, the and (DR)(p^gg) corresponding to the ideal 
system will be referred to as (ref) (ovss) (ref) * 

returning again to 3-12 and 3-13 and solving for the product 

f(A) f'(A), we have for the two equations 



f(A) f'(A) 



2 S(tg) 




Kl 



3-21 



where 



Kl 



^ ^(tg) ^^(cs) 



^6 

f(A) f'(A) = 8 S(tg) S2(cs) I (DR)2( 



OVS8, 



TdrJ^ 



3-22 



ovssj 



where 2 

K - Cd 

2 ■ 8 S(tg) s^(„) I 

We shall now have to assume arbitrary tolerances on the 
undamped angular frequency and on the damping ratio of the system. 
Having made these assumptions, it will now be shown that certain 
definite limits can be set on the non-linear device used in the feed- 
back loop with respect to its departure from an ideal square root 
generator. We shall first discuss the Wn parameter and attempt to 
determine the characteristics of the non-linear device from tolerances 
on Wj^. 

Let W^ of the non-ideal system vary from by some 

fraction of Wy^(j.ef) such as i f ^n(ref)* allowed two 

limiting angular frequencies defined as 

^n(L) " '^n(ref) “ ^ ^n(ref) 
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the lower limiting natural angular frequency and 



'^n(U) ■‘'^n(ref) + ^ '^n(ref) 

the upper limiting natural angular frequency where f is to be defined 
as the tolerance to be set upon the system expressed as a fraction of 
Wn(ref)* 



Thus, we see that for a given tolerance, the right side of 
3-21 is uniquely determined for the two boundary values defined in 
3-23 and 3-24. Solving 3-21 for the slope of the (A, Eq) curve, we 



f'(A) 



Ki Wn^ 

f(A) 



^ ^(tg) 2^(cs) ^(A) 



3-25 



and hence two limiting slopes are determined at every point in a 
plane in which A is the abscissa and f(A) = is the ordinate. We 
define the lower slope at any point (Aj^, f(Ax)) as 



f'(Ai)L 



1 (Wn(ref’) ~ ^ Wn(ref))^ 1 

2 S(tg) [i2;Mj S2(cg) jg.ij 



3-26 



and the upper slope at any point (Aj^, f(Ai)) as 



f'(Ai)u 



I (Wn(ref) ^ 1 

2 S(tg) [i2;M] S2( cs) [e;ij f(Ai) 



3-27 



Fig. 3-3 shows the mapping of the A, f(A) plane for a given 
tolerance of • Kany significant points of interest can be gained 
from a close study of Fig. 3 - 3 * It is first to be noted that the upper 
and lower limiting slopes along a line f(A) = const, are constant. If 
the tolerances on W^^ are increased, that is we are allowed to deviate 
by a greater percentage from throughout the operating range, 

the differences in slopes at any point are increased and hence the 
openings between the lineal elements are increased. Such a diagram as 
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UA) = E. (VOLTS) 




FICURE 3-3 

POSSIBLE CHARACTERISTICS OF NON-LINEAR DEVICE 
IN ORDER TO STAY WITHIN ASSICNED DYNAMIC TOLERANCES 
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shown in Fig. 3-3 is called a lineal element diagramf Paths (l) and 
( 2 ) of Fig. 3-3 follow the steepest slope at any point and hence are 
solutions to the differential equation 3-^7. Path ( 3 ) on the other hand 
follows the lowest slope at every point in the A, f(A) plane starting 
from the same point as (2). Path ( 3 ) is a solution of 3-26. There are 
an infinite number of solutions to each equation, each solution differing 
by a constant. Path (4) shows hov; a path may vary and still give an 
within the given tolerances of the system. Thus, it can be seen that a 
non-linear device which has an input-output characteristic whose slope 
at any point stays within the sp-ecified tolerances of will be a sat- 
isfactory device as far as the dynamics on are concerned. The 
(BR) (Q^gg)tolerances will be discussed in the next section. 

As f(A) — ►-o, f^A) CO and of course it will be impossible 
to obtain a device that will start off with an infinite slope. Thus, 
for low input signals, the dynamic qualities of the system will have to 
be relaxed since the product f(A) f*(A) cannot be maintained within the 
required limits. That is will be decreased and the (^^R)(ovss) 
of necessity, increase. The best we can do is try to obtain as steep 
a slope as possible for low inputs and specify the input below which the 
dynamic tolerances are exceeded. 

We have shown that for the ideal system 

f(A) = Eo = S(sr)[ei;e] si(gg)[A;e] 3-28 

which is the equation of a parabola. Let us re-define the equation of 
the parabola as 

f(A) = Eo “ PA2 3-29 

and denote P(j*ef) constant associated with the ideal system. 

The limits on P are to be determined from the given tolerances on 
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I 



^^^^(0V3S) 



f‘(A) = iPA-i 



and 



f(A) f'(A) = 

Since we have shovm that 






2 S(tg) S^(cs) 
I 



In this case 



r 



y/ 



S(tg) (cs) 



(DR), 






(ovss) * 2 V 2 S(tg) S-^(cs) I f.(A) 



3-30 



3-31 



3-32 



(^^)(oV3s) 2 



I 1 1 

V S(tg) S2(„) I P 



3-33 



Taking differentials in each case. 



Wn = d P 



3-34 



and, dividing 3-34 by 3-32 gives 



d Wn _ ^ 



Wn 



3-35 



also 



■i (0«(OVS») - - S( 



S(tg) S2(cs) I 



d P 

?5" 



3-36 



and finally. 



d (dr) (ovss) - 
(l^«)(ovss 



d P 
P 



3-37 
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From the preceding equations, we have deduced the important 
fact that a fractional change ± f ^x\{ref) the same fractional 

change in the constant associated with the ideal square root curve w^hile 
the same fractional change in (DR)(Qygg), abbreviated 1 f(DR), produces 
an equal and opposite fractional change in the constant associated with 
the ideal square root curve. 

Fig. 3“4 summarizes the results of the previous paragraph. 

Let fi be the specified tolerance on ^^n(ref)> i*®* »^n have any 
value between "^n(ref) ^1 ^n(ref) ^n ” ^1 ‘^n(ref)* Correspondingly, 
the constant P will have any value between F(j.ef) ^l^(ref) ^(ref) ” 
flP(ref)* addition, the (^^^)(ovss) assume any value between 

(ovss) (ref) " ^1 (ovss) (ref ) (ovss) (ref)"^ ^1 (ovss) (ref ) • 

It is important to note here that for a certain given tolerance on the 

natural angular frequency, there automatically corresponds the same 

tolerance on the damping ratio parameter through the upper and lower 

bounds of the square root curves. 

In Fig. 3-4* we have shown how tolerances on the dynamic 

parameters allow an upper and lower bound of square root curves about 

i 

the ideal curve, Eq = ^(ref) remains to connect the development 

in connection with Fig. 3-4 with that of the lineal element diagram of 
Fig. 3-3. We have already shown that 

Eq = f(A) = P d 3-38 

f'(A) = I A-^ 



2 P A2 

- P^ 

■ 2 f(A) 



3-39 
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f(A) = 



^o~<P(ref)+ f |P(ref))A^rot) 




FIGURE 3-4 

LIMITING SQUARE ROOT CURVES FOR A GIVEN TOLERANCE 
ON THE DYNAMICS OF THE SYSTEM 
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Where we have defined P as the upper and lower bounds of 



the square root curves corresponding to the tolerances imposed upon the 
dynamics of the system. That is, P will assume the two limiting values, 



^1 “* ^(ref) ^1 ^(ref) 



3-40 



Thus, since 



Using 3-19 



^2 = ^(ref) " ^(ref) 



^(ref) ~ ^(sr) ^(sg)^ 



P = S(3r) S(3g)S 1 S(sr) S(3g)" 



2 

P = 



Hst) 2 (sg)^ i 2 (sr) S(sg)i 



d2 _ 



W„(r,f) i h «n(ref) 32 ,^^) 



Hence, going back to 3-39, we obtain 

= a ' 7u) 

But, by equation 3-27, we have already shown that 
" 2 S(tg) S(^ 3)2 Vref) ^ Mref) 



3-Ul 

3-42 

3-43 

3-44 

3-45 



3-46 



3-47 



From the fact that 3-46 and 3-27 are identical, we conclude 
that the tolerance limiting square root curves shown in Fig. 3-3 
provide the limiting slopes not only within the band of square root 
curves but automatically give the slopes necessary to map the entire 
plane as shown in Fig, 3*3 • Another important consequence in the 
development of 3-46 is that an assigned tolerance on the natural 
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undamped frequency, simultaneously fixes the same tolerance on the 
damping ratio. 
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CHAPTER IV 



Transient Analysis of the Orthogonal Vector Summing System 



In this chapter, we shall show in detail the response of the 
system to input step and ramp functions. V/e shall take as our starting 
point equation 2-18, and let be unit step function 
Thus, 

^ Eo^ + 2(DR)(ovss) ±_ = 1 4-1 

W7 dt2 ° Wn dt ° ° 

where we have let square of a unit step function applied at 

t = 0. Assume the system at rest at t = 0. 

The response equations to be analyzed iji this chapter will be 
solved by the Laplace transformation method. Tlie Laplace transform of 
a function f(t) will be expressed as 

= F(p) = /' f(t)e“P^dt 4-2 

Taking the Laplace transform of both sides of 4-1, and 
assuming zero initial conditions we obtain 






. aMiovssi p . ij Ep2(p, = 1 4.3 

Here, we have designated E^q(p) ^he Laplace transform of 



Rearranging 4-3, 



[,’2 2 



E^O(P) 



4-4 



p(p2 -f 2(DR)(^ygg) V/^p ^n^) 

There are three cases of equation 4-4 to be considered. The 
three cases to be considered are for (DR)(ovss) ^'^E)(oygg) = 1, and 
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and (^'^)(ovss) ^ 

Before procedin^ further it v:ill be important to note that 
because we are dealing with a second order system with a finite moment 
of inertia, it is impossible by any finite forcing function applied at 
t = 0 to produce a displacement of the shaft which will in turn produce 
an Eo(t). This is of course tx^ue even though we have made the simpli- 
fying assumption that there are no dynamics involved betv/een the shaft 
and Eq. x.ot only will there be no displacement at t = 0 to a forcing 
step function, but the rate of change of response is also zero. Thius, 
the square root of the time resj.onse to be derived from 4-4 will have 
to satisfy these conditions. 

Equation 4-4 when solved by the usual method of partial 
fractions gives as the time response to a unit step function of 

/I - (ovss) / 

“(^^^)(ovss)^ (t^^)(ovss) ^ 

where the djnamic response angle, abbreviated 

^^^>(0V3S) 

2 

To find the initial values of Eq (t) and its first three 
derivatives, we note the following from "Transients in Linear Systems", 
by Gardner and Barnes.^ 

Given a ~ 

Then 

f(0+) = lim p F(p) 

p-*-oc 
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f(0+) = lira 1^ F(p) - f(0+)pj 

P Op 

f(0+) = lim F(p) - f(Os ) ~ ^ r(0+)pj 

P ^ 30 

2 

p^ F(p) - f(0+)p2 - ^ f(0-^)p2 - ^ f(0^)p 

, 00 

2 

Hence, applying the above criteria to E ^(p) developed in 
4-4, we obtain 

Eo^(O) = 0 

jt - ° 

Eo^(O) = «„2 

i ^+2 



3 

f(0+) = lim 

dt^ 



^ Eo^(O) = - 2(DR)( 



ovss) 




We shall abbreviate (^^)(ovss) remainder of 

this chapter to shorten the notation. 

From 4-5, we see that 



Eo(t) = 




1 

\/l - (dr)2 



g-(DR)Wnt sin(W^y^l - (DR)2 t - (DRA) ) 

4-7 J 



Equation 4-7 represents the response of a passive system. 
Thus, even with no damping in the system, the maximum value of the 
second term on the right would be unity. Since a fairly large amount 
of damping will be used for the system, thr 2 second term of 4-5 will 



51 



rapidly become small with respect to unity, and thus enable us to 
approximate EQ(t) closely by two or three terms of a binomial expansion 
of 4-7. 

It is shown in "Differential and Integral Calculus", by 
Granville, Smith, and Longley that the expression 



m 



(1 1 - x)^ = 1 mx ^72 



(m - 1) ^2 



X ♦ 



4-8 



and is valid when m is not an integer providing |xj <1. 

Applying 4-8 to 4-7 we obtain as a good approximation of 

Eo(t), 



Eo(t) = 1 - 



1 - . . 1 3ln(W„ Jl - (DR)2 t - (DHA) 

2 /I - (DR)-^ / 

■ 5 (I - 0>i0' t - (dra)(„,,3)) 



Taking next the case where (DR) = 1, 4-5 becomes 
Taking inverse transforms of 4-10, we obtain 



4-9 



4-10 



^ 2 



(t) = 1 - (w^t . 1) 



Or, finally. 



Eo(t) = 



1 - e“^nt (W^t + 1) 






4-11 



Expanding as before, in a binomial series, 

Eo(t) = 1 - (W^t 1) - (W^t l)^ - . . 

Taking finally, the case where critical damping is exceeded, 
we have (DR) > 1 and hence the characteristic equation associated with 



52 



4-5 factors into 



(p + (DR)Wn - Wny/(DR)^ - l)(p 4 (Dfl)Wn w'ny/(DR)^ - l) 



Thus, 

Eo^(p) = 



Wr, 



p(p -r (DR)Wn - W^_y(DR)2 - l)(p ^ (DR)^^ -r \^(DR)^ - l) 



4-13 



The time expression of 4-13 becomes 
2 {.^ ^ - 1 1 



Eo^(t) = 1 + ^ 



(DR)2 - (DR)y/(DR)2 - 1 - 1 



[-(DR)Wn + Wn^/(DR)^1 ' 



+ i 



(DR)2 + (DR) _^(DR)2 _ 1 - 1 
Thus we have finally, for (DR) > 1 



^[-(DR)W„ - Wn^(DR)2 - 



t 

4-14 



Eo(t) = 



1 + k 



+ k 



(DR) 2 - (DR)-^(DR)2 - 1 - 1 



(DR)2 + (DR) ^ (DR)2 _ 1 _ 1 



[-(DR)W^ ^ Wj,^(DR) 2 - 1_ t 
[-(DR)U„ - v/„^(DR) 2 . 1 ]t 



(DR) > 1 4-15 

It is not feasible to expand this expression by the binomial 
approximation . 

Curve I of Fig. 4-1 shows the solution of 4-5 for (DR) = .1. 
If the square root of the values representing curve I are taken at 
every point, we have the exact solution shown as curve II, and which 
represents 4-7* Since this latter result is difficult to obtain 
analytically, we have expanded 4-7 through the binomial series approxi- 
mation to three terms in 4-9* Plot III sho;>rs the approximate value of 
Eo(t) to a unit step input function using only the first two terms of 
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FIGURE 4-1 

NORMALIZED TRANSIENT RESPONSE CURVES FOR 



4-9* It/ is to be noted that the approximation is very good except 

for the initial period of ^ = .1. Since damping ratios higher than 

^n 

.1 are actually to be used, the approximation will be even closer to 
the actual solution except that as the damping ratio increases the 
time for the initial overshoot will correspondingly lengthen and hence 
the initial large deviation will last for a longer value of 

Throughout this thesis, it has been emphasized that the 
terminology used such as damping ratio, and natural frequency refers 
only to the linear second order system relating the squares of the 
output and input functions. However, one must be careful and not 
apply this terminology to the specific response associated with the 
input and output functions themselves. 

In this presentation, we shall use the logarithmic decrement 
as our basic description parameter for the system for (DR) < 1. Loga- 
rithmic decrement will be defined as the natural logarithm of the ratio 
of the amplitudes of two values of the response to a step function 
differing by a period as for example two successive maxima of oscil- 
lation as shown in Fig. 4-1. 

Referring to Fig, 4-1, we let and A 2 be the amplitudes of 
two successive peaks of the linear second order response of the squared 
functions while and A 2 * are two successive peaks of the response 
for signal input as a function of signal output. 

Hence, the logarithmic decrement for the linear second order 

A. 

response, abbreviated (LD)(lin) “ 4-l6 

From Equation 4-5, we see that the instantaneous height of 
the overshoots shown of Fig, 4-1 is 



A . (ni)i sin(W„^l - (DR)2 t - 



) 
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and hence using 4-16, 



• -(DR)Wnti 

Umj _(dR)W t2 



4-17 



since and A2 were taken at the peaks and ^ = t2 - t^j^, the period of 
the damped oscillation. 

From A-I7, we obtain 

(LD)(ii^) = (DR)W^ T where T is the time for two positive 

, thus 

1 - ^UR;‘- 
"V 



,,.m_27r_ 2 7T 

successive peaks. But T - -rp “ ~ — / ' , s' 

W Wn (DR)' 



_ + 2 7T(DR) 

(LD) lin) - / ^ . 

VI - (DR)‘ 



4-18 



We shall next define the logarithmic decrement for the actual 
response of the system to ^ step input as, 

4-19 



<'-“)(ovss) = =ilogii 



*2 *2 
From 4-19, we have immediately that 



(LD)(ii„) = 2(L0)(,„,3) 4-20 

Using 4-20 in 4-17, we obtain the important result that 

(i-W(ovss)= '•-21 

yjl - ( DR )2 

This last result, gives a means of describing the non-linear 
response in terms of a parameter describing the linear squared system. 

Up to now, the development of this chapter has been with the 
ideal system which assumes a perfect square root function in the closed 
loop. The last chapter showed how close a square root function had to 
be generated in the closed loop in order to stay within certain pre- 
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assigned tolerances in the damping ratio and the undamped natural 
frequency. Returning to this line of argument, if we have a change 
in the damping ratio, d(DR)(Q^gg) then from U-20 we obtain 

d(LD)( ) = ^(1 (DR)2) (DR)< 1 4-22 

Applications of the response of the system to a forcing step 
function will be made in connection with the laboratoi^'* test system. 

The constancy of the dynamics for various magnitudes of the input signal 
will be tested by applying small step function increments upon the dif- 
ferent static operating levels and recording the output response. If 
we have a perfect square root generator in the feedback loop, the 
response to a step function increment will not depend upon the magnitude 
of the forcing function, and hence a non-dimensional response curve can 

be shown by plotting — ^2^ for any point of the operating region. How- 

^in 

ever, we shall use an imperfect square root generator in the feedback 
loop and hence small increment forcing step functions will be used and 
the operating points specified. Chapter VI will show in non-dimensional 
form, such responses as obtained from a Sanborn recorder. 

The previous analysis with logarithmic decrements was made on 
the basis that the (DR)(QYggj <1. We have not developed an easy para- 
meter to describe the non-linear system in terms of the linear squared 
system for {^^){ovss) ^ However, the response can be obtained 

fairly easily from 4-12 or 4-15* 

It was felt that to make this chapter on transient analysis 
complete, an analysis of the system response to a ramp forcing function 
should be included. Much study was devoted to the comparison of the 
response of a regular second order system to the one presented here, but 
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the results were disappointing because no simple relation or design 
parameter could be found relating the two. However, a brief develop- 
ment of ramp responses follows for (DR) < 1. 

Taking the Laplace transformation of 2-lS, and assuming zero 
initial conditions, we obtain 



Wn 



or 



Eo'^(p) = 



Ein^p) 



l,p2, ^P^l 



W, 



4-23 



In developing the response to ramp functions, it will be 
best to work from the point of view of errors (both steady state and 
transient) between the input and output functions. 

Treating 4-23 first as a linear second order system in E 

we define the error for such a system as Cl |e\„(p) - E?„(p) 
where 



Eq^(p) - E^in(p) 





-f 



2(Mp + 1 



or 

Eo^(p) - E^in(p) 



^\n(p) [p^ ^ 3(DR)’^'n P] 
p2 -h 2(DR)Wp p 



If (t) 
in' ' 



a unit ramp function = t, 



4-24 



4-25 
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and the linear second order error equation becomes for DR < 1 

P P P 2(DR)Wn 

Eq^(p) " %n^(p) = ~7~2 3 ^ 

p(p^ + 2(DR)Wn p 

p 2(DR)Wn 

~ p(p (DR)Wn - jWn^l - (DR)2)(p ^ (DR)Wn + jWn^l - (DR)^) 

4-26 

- ^ ^3 ^ 

P p (DR)Wn - JWn^l - P (DR)Wn 0W„^1 - (DR) 2 



where 



Ko = 



Ki = 



(p) (Eo^ - Ein^) 



_ 2 (DR) 



W. 



n 



p = 0 



(p (DR) - jWn^l - (DR)2 )(Eo2 - E^^^) 



. . -1 2VT^(DR)^ 

»J tan Q5PQ 



p = - (DR)Wn - (DR)' 



2 JWn^l - (DR)' 



Hence, the complete time expression for the linear second order 
error response to a unit ramp is 

* Wn^l . (DR)i * * * " ^ ■' 

0 



(DRA) 



(ovss) 



4-27 



where 



(°'^)(ovss) = 



,-i gV^r^DRjl. (DR)< 1 
(DR) 
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then 



If, on the other hand, we let £• - a unit ramp function = t. 






2 



and 4-25 becomes 



Eo^(p) 



Ein^(p) 



2 I 


[p * 2(DR)W^] 


1 


p2 


p^ 


2(DR)WnP ■* 





2I 


[p + 2(DR)Wn. 


1 


p2 


P (DR)W^ - - (DR)2j 


[p (DR)Wn ^ jWny 


- (DR)^ 



(DR) < 1 4-28 

Expanding by the same method of partial fractions as in the 
previous case gives 






Wn2(2(DR)2 . 



(DR)' 



■(DR)Wnt sin(Wj,yCr^R)2 t ^ (DRA)^^^^^)) DR < 1 



4-29 



where 



(DRA)/ n = tan' 
(ovss) 






An interesting result from 4-29 is the fact that if the 
(DR)(ovss) is made equal to .5, then the first steady state term in 
4*29 is eliminated* However, the second term shows that the error 
between the squares of the input and output builds up linearly with 
time. 
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CHAPTER V 



Steady State Harmonic Analysis of the 
Orthogonal Vector Summing System 



An important method used to determine the dynamic charac- 
teristics of a general electro-mechanical system is the method of 
steady state sinusoidal harmonic analysis* Amplitude and phase graphs 
are plotted as functions of forcing sine wave frequencies* In this 
method, all transient responses are considered as having died out 
leaving only the steady state conditions. The dynamic quality of the 
system is then determined by the pre-assigned limits set upon the 
magnitude change and the phase change in comparison with some reference 
magnitude and phase determined at some reference frequency, the latter 
usually being zero. 

In the system being considered here, it is impossible to 
follow such a procedure because of the squaring properties of the 
torque generators. Thus, only the absolute value of a function acts 
on the system, and of course |cos Wt|, which looks like the following 
function, 




has too abrupt a change to be of any value in evaluating the speed of 
response of the system. 

If, on the other hand, we impress upon a constant input 
a function of Wt, Fq_ cos Wt, where Fl = Fo> the discontinuity can be 
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prevented. Thus, let Bin = F© ♦ ^os Wt. Since the equation giving 
2 2 

Eq in terms of is linear, all of the solutions for inputs of E^ 

will be found through the squared functions. The equation for Eq^ in 
2 

terms of Ej[n is 



1 d^ 
Wt? ^ 



at 



W, 



n 



5-1 



An important significance of this linearity is that if E^ 
consists of a sum of terms, a solution of the equation can be obtained 
for each of the terms independently and the overall solution obtained by 
superposition of the separate solutions. 

Taking an input function, Ej^j^ = Fq cos Wt, and squaring 



gives 



Ein^ « Fo^ 2 Fo Fi cos Wt f Fi^ cos^t 5-2 

From trigonometry, we have the following relationship: 



cos^Wt = i i cos 2 Wt 



5-3 



Thus, 



or finally. 



= Fq^ ♦ 2 Fq Fj^ cos Wt 



« 2 p ^ 

El n 



cos 2 Wt 



Ein^ = Fq^ + - — + 2 Fq F^ cos Wt + - — cos 2 Wt 5-4 
2 ^ 1 ^ 

Taking the Fq + ^ — term and inserting it into equation 5-1 

2 Ei2 

gives the response Eq = ^ the steady state. It is to be 

stressed that this harmonic analysis deals only with steady state 
responses of the system. 

The steady state solution of equation 5-1 to the fundamental 
sinusoidal forcing function will now be shown. 
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Rewriting equation 5-1 as follows: 



1 ^ Eo^ . 2iM 



Wn^ dt^ 



W, 



'n 



^ Eo^ -h Eo^ = 2 Fo Fi cos 
dt 



- Re 2 Fq ( 



Let 



z = Eo'" 

Assuming a particular solution of the form, 

Zp = A ejWt 

Zp = jWA 

Zp = - W2 A ej'* 

Upon substituting the above functions into equation 5- 



we obtain 



- id- ^ 1 + j 2(DR) iL 



W, 






A = 2 Fo Fi 



A = 



2 Fq Ft 



ildy ^ 1 -r j 2(DR)W- 









Thus 



Eo = Re 



2 Fp 



1 - w2 + j 2(DR) iL 

W„ 



w, 



Eq^s Re 2 Fq Fi(co 3 Wt -» j sin Wt) 
1 - (FR)2 + j 2(DR)(FR) 



where 



(FR) = ^ 
'^n 



Wt 5-5 
5-6 



5-7 

• 6 , 



5-8 

5-9 



63 



_ 2 Fq Fi(cos Wt j sin Wt)(l - (FR)^ - j 2(DR)(FR)) 

[l - (FR)2j ^ 4 [2(DR)(FR)j ^ 



or 



LEo£i , 

jl-(FR)2 . r2(DR)(FH)J' 

or finally 

^2 2 Fq Fi 

1-12 



- (FR)2J2 ^ |"2(DR)(FR)‘j 



(1 - (FR)^)cos Wt 4 2(DR)(FR) sin Wt 

5-10 

cos (Wt 4 (DRA)i) 5-11 



where 



(OR*), = 



(FR)' 

On a similar manner, the response of equation 5-1 to a 
Fi2 

function cos 2 Wt is of the form 



Fl^/2 



v/[l - (FR)i^] 2 4 [2(DR)(FR)i]‘ 



cos (W^t 4 (DRA)2> 5-12 



where 



(DRA)2 = tan~^ - 



Sl.ce (FR)i= S 



2(DR)(FR)i 
1 - (FR)i^ 

2(FR) 



Thus, 5-12 becomes 



F1V2 



■\/[l - 4(FR)2]2 + |^4 (br)(FR)] 



cos (2Vft 4 (DRA)2) 5-13 
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where 



(DBA), = ta„-l - idmsm 
1 - /.(FR)^ 



In the case of 

4(FR)2« 1 

(DRA)2 = 2(DRA)j^ 5-14 

This can be seen from the fact that for (FR) « 1 the tangent 
of the angle can be replaced by the angle and (DRA )2 — - 4(DR)(FR) as 
compared to (DRA)i^ = - 2(DR)(FR), 

This approximation is good even for (FR) - .3 • To show this 
fact, we shall let (DR)(ovss) “ ^ 

(DRA)i = tan“l — 

1 - .09 

= - 33.4° 

Letting (DRA )2 = 2(DRA)j^, gives - 66,8° for (DRA )2 
Actually, 

(m)2 = ta„-l = -62« 



Thus, corresponding to 



E. 2 - F 2 ^ F 
*in - ^o +2 



♦ 2 Fq F]^ cos Wt + — cos 



2Wt 



5-15 



we have in the steady state 



2 - IT. 2 



Er>^ = F, 



Fi< 



0*2 



2 Fq Fi 



\/ 1 - (FR)2]2 4 [2(DR)(FR)]2 



cos (wt -f (DRA)j^) 5-16 



Fi^/2 



\/|l - 4(FR)2] 2 [4(DR)FR)] 



cos (2Wt (DRA) 2 ) 



65 



where (DRA)]^ and (DRA )2 have been defined in 5-11 and 5-13 • 

It is now necessary to relate in some fashion, the squares 
of the inputs and outputs to the actual inputs and outputs themselves. 

Starting with the expression Eq 2 - Ein2 = (Eo - E^^) , 

we can make the approximation 

Eo^ - Ein^ = 2 Ein(Eo - Ei^) if Eq « Ein- 5-17 



This approximation will be valid because we shall restrict 
ourselves to small deviations of Eq. 

2 

Dividing both sides of 5-15 by gives the result 




Hn 




5-18 



We shall define the deviation of the performance ratio of 
the orthogonal vector summing system as 



D(PR) 



(ovss) 



4n 



- 1 = ^ 




5-19 



Eq^ 

where 5-15 and 5-16 are to be used to give the expression for r — o - 1. 

^in 

Information from 5-19 will be broken into two parts, each 
part conforming to certain specifications. In this particular treatment 
we shall use the damping ratio associated with the linear second order 
performance equation relating the squares of the input and output 
functions as the necessary parameter to give us a family of performance 
ratio curves plotted as a function of frequency. The first of these 
aforementioned parts will be a comparison of the ratio of the absolute 
magnitudes of and while the second part will give the phase 
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lag of Eq as compared to a functional forcing frequency. It 

is to be noted that the amplitude deviation will be a comparison of 
the maximum amplitudes of input and output quantities. From the infor- 
mation obtained from the ratio of the absolute magnitudes of and 
Ein^ we can then use 5-19 to find D(PR)^Q^ggj as a function of frequency. 
In the expression for that is E^^ = ^o ^1 

and F-^ have been specified as not necessarily equal. This general 
approach has been carried through with an eye on the future applications 
of this type of analysis. In considering the quality of response of 
the system, Fq will be made equal to Fi* However, in considering the 
effect of noise on the system, F]_ will usually be a small fraction of 
Fo* We shall apply our analysis only to measure the quality of response 
of the system, considering a study of noise beyond the scope of this 
thesis. 

Thus, letting 

Ein = Eq -► Fq cos Wt 5-20 

Then 

^in = 1 + cos Wt 5-21 

Fo 

Let Fq be redefined as Ein(ref) 

- = 1 + cos Wt 5-22 

%n(ref) 
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Then, 



■•in 



%n^(ref) 



= 1.5 ■» 2 cos Wt + .5 cos 2 Wt 



5-23 



After substituting this quantity into equation 5-1, 



= 1.5 t 



E in(ref) ^ , [2(DR) (^^ 33 ) (FR) ] ^ 



cos(Wt (DRA)]^) 



V[i - ).(fb)=]2 , [udr)(„,,^)(fr)] 



cos( 2 Wt ^ (DRA) 2 ) 



5-24 



where, for the sake of completeness 



(DRA), = tan“^ - 



4(DR)(FR)^ 

1 - 4(FR)‘ 



Fig. 5-1 shows how the input function ■ ^ of 5-23 is 

k 2 . ^ in(ref) 

compounded. It is to be noted that has a frequency, W, but 

^in(ref) 

the required addition of the second harmonic causes distortion as shown 



in Fig. 5-1* For sufficiently high forcing frequency, — 

17-2 E^in(ref) 



will 



lag the input as shown and the absolute value of 



e2 . 



will drop 



'in(ref) 

off in general or it could very well peak for low damping ratios near 

the natural frequency of the system. 

Fig. 5-2 shows in non-dimensional form a plot of the 

D(PR)(ovss) a function of forcing frequency using equation 5 - 19 * 

e 2 • 

The absolute value of can be seen from equation 5-23 to be 4* 

E2in(ref) 
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FIGURE 5-1 

COMPOSITION OF THE NORMALIZED SQUARED 
IMPUT AND COMPARISON OF ^ TO 
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lao 



can 



can be calculated from 5-24 with (DR)(ovss)^® parameter. 

^ In(ref) 

The second harmonic component will cause a resonant point at (FR) .5 
in addition to the usual resonant point, hence the broadening of this 
peak over the usual form in which it is seen. The fact that for low 
damping values the resonant peaks become narrow and high, brings out 
an interesting point in connection with the curve for = .1, 

Here we have two distinct peaks, the first due to the second harmonic, 
coming at a frequency ratio slightly less than o5« 

In order to find the phase response of the system to the 



input function 



-•in 



= 1 cos Wt, it will be best to look at the 



Ein(ref) 

squares of the outputs and inputs in rotating vector form. In setting 
up the relationship between the input and output rotating vectors, let 
us assume (DRA)]^ = 30^ and (DRA )2 = 60^ as a typical case. Further, 
let us assiime for simplicity that the amplitudes of the response vectors 
are the same as the corresponding input quantities. 

It is obvious that the phase angle between input and output 
will not remain constant throughout a complete rotation of the funda- 
mental vector due to the second harmonic vector rotating twice as fast. 
This changing of phase angle throughout one period becomes more pro- 
nounced as (DRA )2 builds up* The best that we can do for the given 
system is to predict where the (DRA) starts to build up to a few 
degrees. Inasmuch as the fundamental vector is four times the second 
harmonic for most cases (not near resonance for low damping ratios), we 
can, for small phase angles, approximate the total phase from the funda- 
mental as can be seen from part I of Fig. 5-3* also must define this 
phase angle to be determined at t = 0, 2 , etc. as shown in part I of 

5-3 • f'ig* 5-4 gives the phase lag of the fundamental as a function 
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= CONSTANT INPUT-OUTPUT 



assume: 



= input (DRA)| = 30® 

* output (DRA)2= 60® 




4 




INSTANTANEOUS RELATIVE POSITIONS OF^ OUTPUT AND INPUT FUNCTIONS 
FOR THE FORCING FUNCTION. ||^=i.5^-2coswt+.5cos2wt 

in(ref) 
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100 



FIGURE 5-4 



(DRA)(ovss) AS A FUNCTION OF FREQUENCY FOR 
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of frequency. Caution must be taken to use Fig. 5-4 only for small 
angles. 

In concluding this chapter on the steady state harmonic 
analysis, it must be pointed out that approximations have to be made 
due to the non-linearities of the system. These approximations give 
good results at the point where the magnitude of the (D)(PR)(^ygg) 
just start to build up and where the phase lags are only a few degrees. 
We have assumed in order to use the material presented here that a 
non-;linear device has been inserted into the feedback loop so as to 
keep the dynamics essentially constant. 

This chapter emphasizes the difficulties in obtaining a 
method of analyzing a non-linear system so as to give a clear cut 
picture over the entire operating range. 

It is felt by the author that a step function response type 
of analysis is the more suitable method of ascertaining the speed of 
response for this type of system. The experimental results presented 
in the next chapter will bear out this opinion. This chapter was 
included mainly in keeping with the underlying spirit of this thesis, 
namely that it be a study of the dynamics associated with the given 
system. 
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CHAPTER VI 



Design and Results of a Laboratory Test System 



The general theory developed in Chapters II, III, and IV 
will now be applied to the design and testing of an actual system. 
Because of the fact that exact specifications of microsyn units are 
still considered restricted, we Shall not use the actual values, but 
we have chosen the specifications so that the test results do actually 
give a true picture of the operation of the system. 

Thus, using Fig, 2-2 as our starting point, we shall use the 
following design parameters for the system to be designed in this 
chapter, 

S(tg)[i2;M]= 6-1 



S(sg) [A;i,n,e] = 3-25 x 10'^ 



volt 

rad ma cycle/sec 



i(ex) = 

^(ex) “ cycles/sec 

^(rot)(max) " radians 

l(tsm) = 330 gm cm2 

Using the above excitation values for the signal generator, we obtain 

volts 



0(sg) [Aje] 



= 13 



rad 



6-2 



Since 50 ma is a desirable upper limit of driving current 
for each torque generator, the maximum torque produced per torque 
generator is 
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M(tg)(max) = ^(tg) ^ (max) 

^ 1^2 M X 2500 ma^ 

ma^ 



3 106, COO dyne cm 6-3 

Hence, the total maximum torque produced by the maxium input 
signals will be 212,000 dyne cm. Thus, the maximum torque produced by 
the feedback voltage Eq will have to be 212,000 dyne cm. 

Fig. 6-1 shows in schematir' form the torque generator winding 
arrangement, i^ and iy are connected in single axis arrangement on 
separate microsyn units, i^, the current due to the feedback voltage 
Eo, is connected in single axis arrangement but windings of poles 1 and 
3 of each microsyn unit are connected in series so that the feedback 
current io produces twice the torque that it would if it were to actuate 
only one microsyn unit. This is necessary because we need to keep the 
maximum torque generator current at 50 ma. 

Referring again to Fig. 2 - 2 , it can be seen that in the steady 



state. 



(cs)]^[ejiJ ^(tg)j^ ^ ^ (cs) 2 [e;ij '^(tg )2 ^ ^ (cs) (tg)^ 

6-4 



. M2 = 

.2 



i.’ 2 



Since 

and 



^(cs)i ~ ^(cs)2 
^(tg)l “ ^(tg)2 



S2 

(cs)i,2 



'(tg)l ,2 




«2 q p* ^ 

( 03)3 (tg)q O 



6-5 



Hence o o o 

Eo'" z Ex'" if 
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I 
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WIRING ARRANGEMENT OF THE TORQUE GENERATORS USED IN LABORATORY SYSTEM 



6-6 



[.®5iJ ^(tg)l,2 “ ^ (cs)3 [eji] ^(tg)3 

For the system described here, we can define a fictitious 
third torc^ue generator as shovm in Fig* 2-2 which will have twice the 
sensitivity of a single torque generator. The resistance of this series 
connection of torque generator windings will be double that of single 
microsyn operation. However, because of the fact that a current source 
actuates the torque generators, a change in resistance in the windings 
does not change the driving current since current sources have infinite 
internal resistance. 

Referring to 6-6 and noting that 

^(tg)3 = 2 , 

^ (cs) 3 [e;i] ~ ^ ^ (cs)i [e;i] 

S(cs)3 ~ ^(cs)i [eji] 

From equation 2-8 we have defined 

S(ovss)[A;M] = 2(tg) [i2;k] ^^(cs)[eji] ^^(sr)[e2;e] ^(sg)[Aje] 

where it has been assumed in 2-8 that all the torque generator and 
current source sensitivities are respectively equal. However, in the 
laboratory system we have defined |^j_2 j = 2 2 [^^5^] 

Nevertheless 2-8 is still valid here since we have kept the product 
S(tg)[i2;K] S2(cs) [e;i] constant. 

It is desired to use the full 10° or .174 radians for rotor 
angle in designing the system, this angle being produced at the maximum 
torque level. 

Hence 
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^(ovss) 



212,000 dyne cm 
•174 radians 



= 1,220,000 dyne cm/rad 

Letting Ex(jnax) > -^(max) 50 volts each, 

=(«)X = - 1-Aoit 

^(max) “ ^(max) “ 50V 

Eo(max) = J 25 OO * 2500 = ?0*7 volts 



Hence 

3 _ i 0 _ ma( _ m^) . ma/volt 

70.7 volts 



This latter value checks with the result of 6-7- 
Using equation 2-8, we obtain 



S2 



(sr) 




^(ovas), [a;H] 

^(tg)3j^i2;Mj ^^(03)3 j^e;i] ^(sg) j^Ajej 



S2 



(sr) [s^jej 



dyne -cm. 



1 , 220,000 



.848 X 102 x .5 ^ x 1.3 ^ x lol 



iha"= 



volts 



rad 



_ 2210 

volt 



or 

2(sr) j^e^;ej 



47.0 



volts 

volt^ 



Figure 3-1 shows the characteristics of the square root generator 



6“8 



6-9 



6-10 



6-11 



6-12 



6-13 
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To actually obtain a device that will have a performance 
characteristic similar to Fig. 3-1 is very difficult. This is especially 
true near the origin where the slope of the curve becomes infinite. 
However, through the use of vacuum or selenium rectifier switching cir- 
cuits in cascade with a linear amplifier, the curve of Fig. 3-1 can be 
approximated to a high degree except at the origin. The latter part of 
the chapter will be devoted to the design of a suitable attenuator which 
will stay within certain specified tolerances of the dynamics of the 
system. 

Fig. 6-2 shows an electrical schematic of the system from 
which the experimental results will be obtained. High gain DC ampli- 
fiers with unity feedback, feeding into parallel connected 6SIf7 triodes 
in each case act as current sources. 

The DC amplifier arrangement as shown in Fig. 6-2 can be 
considered current sources because of the following conditions. The 
gain of the DC amplifier is very high with a phase inversion and the 
input impedance is also very high. Thus referring to the figure below, 
the input to the DC amplifier (point G below) is driven very close to 
ground because of the very high negative feedback. 




Furthermore, since the input impedance to the amplifier is very high, the 
current through Ri equals the current through R2. But since Cl is ideally 
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FIGURE 6-2 

ELECTRICAL CIRCUIT DIAGRAM OF LABORATORY TEST SYSTEM 



at ground, we can consider the current through as depending only 
on Ejj^ and Ri leaving R 2 to vary at will. However, through current 
considerations 



Ri 




or 



Eq _ ^ 

Ein h 



6-14 




actual system we have fed Eq into a cathode follower and the voltage 
at the feedback point of the cathode is still determined only by the 
ratio of the feedback resistor to the input resistor. The cathode 
resistor is then calculated to give a certain desired current for the 
given input voltage. Of course, we have to assume that the output 
tube operates below any non-linearities due to saturation. The cathode 
resistors are determined so that maximum currents of 50 ma each are 
delivered for each respective maximum voltage input. 



about 2.25 volts nns (400 cps), an AC amplifier has been provided with 
a gain of around 500. Tliis amplified signal generator voltage is then 
full-wave rectified and filtered before being fed through the non-linear 
device as DC. The output of the non-linear device gives the desired 
solution, Sq. The components within the dashed line box provide the 
S(sr)|^ei-ej shown in Fig. 2-2. Because of the fact that 

the maximum forcing frequencies will be on the order of a few cycles, 
the square root section will be assumed to add no additional dynamics 
to the system. A VR-75 regulator tube is used aroun*^ the feedback DC 



Since the maximum signal from the signal generator is only 
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an^lifier in order to protect the amplifier from saturation when testing 
the system in the unstable region described in Chapter II. 

We shall now apply the development of Chapter III to the actual 
system just described taking into careful consideration the need for the 
stop near the null in order to maintain a stable operating range. 

All of the operating components of the laboratory system have 
been specified with the exception that the square root generator was 
left ideal and specified by the equation, 

Eo = 47.0 E(gg)5 

In Chapter III, we have denoted the ideal case as 

Eo = f(A) = P(pgf) A5 6-15 



where 

P(ref) 

Numerically 

^(ref) 




47.0 



volt 
— 77T 

V0lt2 




volts^ 

rad2 



or 



volts 

"(ref) - 



6-16 



In setting up a workable system, let us arbitrarily set a 
tolerance on the dynamics so that the natural undamped frequency and 
the damping ratio do not vary from '■^n(ref) and (DR) (qvss) (ref ) 
than 25!^. This tolerance may seem somewhat large, but it will be seen 
that the principles underlying the theory developed in this chapter are 
more clearly and easily shown if a fairly large tolerance is assumed* 

The expression for the natural angular frequency of the ideal 
system as given in 2-15 is 
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- ^^(ovss) A;m] 



where j = 1,220, OCX) dyne cm/rad and I = gm cm*^. Using these 

values 

^n(ref) “ rad/sec 6-17 

Thus referring to 3-23 and 3-24, we obtain 
Wn(L) = 60o8 rad/sec - ( * 25 ) (60.8) rad/sec 

Wn(L) = 45.6 rad/sec 

and 



6-18 



Wn(u) = 60 o 8 rad/sec ^ (« 25 )( 60 , 8 ) rad/sec 

Wn(u) = 76.0 rad/sec 6-19 

The two limiting slope equations are thus obtained from 
3-25 and will be defined as 



f*(A) = 3.91 volts sec 



2 2 ^n 



f(A) 



6-20 



Hence, for the case in question, the upper and lower limiting 
slopes become 

"''“'I, 



r'(A), = 



‘ '*■'11 f(A) 



From 6 -I 5 and 6-16, we obtained the ideal square root 
equation as 

volts 



6-21 

6-22 



Eo - f(A) = (170 ) Ai 



6-23 



Using the arbitrary 25% tolerances, we obtain the limiting 
square root curves as shown in Fig. 3-3 as 
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f(A)L 

or 

f(A)L 

similarly 

f(A)u 



170 - (.25) 170 



volts 

radi 



Ai 



127.5 



volts 

rads 






212.5 



volts 
1 A* 
rads 



6-24 



6-25 



Taking the functional derivatives of 6-24 and 6-25 and forming 
the respective products f(A) f'(A) we obtain 

6-26 



^ f(A) 



a 



22.600 

f(A) 



6-27 



Equations 6-26 and 6-27 check with 6-21 and 6-22 and hence are valid 
for obtaining the limiting slopes throughout the A, f(A)‘ plane. We 
shall abbreviate \y;ot) A to simplify the notation. 

Fig. 6-3 shows in detail the limiting curves of the square 
root generator with the assumed tolerances of 1 25^ using the full 10^ 
of shaft angle. The limiting slopes are determined at any point A, 
f(A) by f(A)j^ and f(A)^. However instead of fixing the curves at the 
origin, let us fix the upper end point at the ideal values corresponding 
to radians which gives an ideal Eq = 70.7 volts. 

If we start at A(m^), Eo(max) always follow the lowest limiting 
slope at every point, we notice that it is possible to stay within the 
dynamic tolerances and Intercept the Eq axis at a value of 47 volts. 

This path is of course trivial because it would be impossible to form 
the right triangle addition for any value of the hypotenuse less than 
approximately 47 volts. Starting again at ^o(max) this 

time following the highest limiting slope at each point we find a much 
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more useful situation. In this case, it is to be noted that it is 
possible to compute Eq down to the limiting value of zero. We observe 
in our case that for A » ^064 rad., Eq = 0. The fact that it is possible 
to allow a rotation of the rotor before allowing an output suggests the 
feasibility of using an electrical stop at a small positive angle in 
order to overcome the condition of instability at the origin described 
in Chapter II. 

It is to be emphasized again that the loop gain does not affect 
the steady state conditions of the system except to change the angle of 
the shaft, but does effect the dynamic qualities. Thus, doubling the 
loop gain from that amount shown in Fig. 6-3 reduces the maximum angle 
to 5^, increases the natural undamped frequency by ^ 2 , and decreases 
the damping ratio by ,^2. 

The first condition can be seen from the following considerations. 
From Chapter II, we have shown that one version of the performance 
equation of the system is 

^(rot) A + Cjj A ♦ 2^(cs) ^ (sr)^e2;e] ^(sg) [A;e] ^ 

= ^ (cs)[e;^ ^(tg) ^ in ^ 



Thus, in the steady state, 



^(max) 



^ in(max) 



/ spl \P "I 

(sr) [e2;ej (sg) |^A;ej 



6-29 



Since the square root sensitivity v;ill consist of a linear 
a-c amplifier demodulator combination in cascade with a passive attenu- 
ating network which will, in the ideal case, follow a square root function. 
Hence, we may vc'rite 
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(sr)[e 2 ;e] ” ^ (amp)[e;e] ^ (att)[^e^;e] 



and 



^(max) ~ 






(amp)[^e;ej ^ (att)[e2;e] ^(sg)[A;e] 



6-30 
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But, since 



2 2 2 1 c 

^(ovss)[a;m] ~ ^(tg)[i^;M] ^ (cs)[e;i] ^ (amp)[^e;ej ^ (att)j^e^;aj (sg) [Ajej 

6-32 

. _ ^ in(max) ^(tg)[j^;M] ^ (cs) [e;i] 

^(max) “ 

O(ovss) I A;M J 

or finally, 

const. 
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A (max) = 



’(ovss) [A;Mj 



6-34 



Thus, we have the two results that increasing the sensitivity 
of the amplifier-demodulator in cascade with the non-linear attenuator 
by the factor yj 2, doubles the loop gain and cuts the shaft angle in 
half. 

Preferring again to Fig, 6-3, we have shown that at any point 
in the (A, Eq) plane 



^ . 8130 

dA L ■ f(A) 



6-35 



and 



dA 



22,600 

f(A) 



6-36 



in order to stay within our 25$ tolerances on the natural frequency and 
dcunping ratio of the system. 
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starting at the point, (lO®, 70.7V), 6-85is represented by 
curve 4 and 6-J56 is represented by curve 3. It will be easier to design 
a non-linear device if we plot 6-36 and 6-36 as shown in Fig, 6-4. As 
long as we stay within the curves (A) and (B) as Eq is changed, then 
we are meeting the required tolerances on Wj^ and (l^R)(ovss)* 
laboratory test system, we shall use a straight line attenuator which 
will abruptly bend at points where the slope of the original line will 
cause the tolerances on the dynamics to be exceeded if carried any 
further. 

Thus, starting at the point, (70.7V, 320V/rad), on Fig. 6-4, 
we proceed along the constant slope III until is lowered to 25.5V 
where we must break in order to stay within the required boundary 
curves A and B. The same procedure is carried out for slopes II and 
I. Slope I is carried straight through to the Eq = 0 line and the 
dynamic tolerances are exceeded at Eq = 3*2V. Because of the infinite 
slope required at the origin, the best we can do is to make the slope 
for low signal levels as steep as possible and specify for what levels 
the tolerances on and (ovss) exceeded. Fig. 6-3 shows how 
curves I, II, and III look on the A, f(A) plane. An important point 
to be brought out here is that slope I crosses the abscissa at some 
small positive rotation of A. Thus, we have achieved our goal of 
putting a stop (electrical in this case) so as to keep A from passing 
into the unstable region. Our next step is to instrumentize a suitable 
system to give the response represented by slopes I, II, and III. 

Our first step in developing an actual circuit that will 
represent the three different regions is to determine from Fig. 6-3 
the actual gains that are required for each region. What will be done 
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OOO'OI 




o o 

o - 



(V),i = 



VP 

®3P 
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Eo = f(A) 



will be to attenuate the amplifier gain required for step I when 
passing into regions II and III. 

Taking the slopes from Fig* 6-3 which will give us the full 
ten degree of shaft rotation for maximum Eq we have 

Slope ^(amp) E;e 

I 192 o3 v/v 

II 68*1 v/v 

III 23.9 v/v 

If we desire five degrees for maximum shaft angle, then each 
of the above gains are multiplied by the square root of two. 

The experimental results to follow will also show how the 
dynamics vary quite radically when a purely linear device is used in 
the feedback loop. For ten degree operation this characteristic is 
represented by slope 6 of Fig. 6-3 and C of Fig* 6-4. Frdm Fig* 6-4, 
we see that the required 25 ^ tolerances are met for values of Eq 
ranging from 20V to 56V* For values of E^ greater than 56V, the 
slope of the linear feedback device is too high as can be seen from 
Fig. 6-4 and hence the natural undamped frequency is more than 25^ 
greater than and conversely the damping ratio becomes less than 

(^^)(ovss)(ref) " ^5^ (°^4vss) (ref ) * opposite is true for 

0 = Eq = 20V. 

Fig* 6-5 shows the attenuator used in the laboratory test 
system which, in conjunction with a conventional a-c amplifier and 
demodulator, will follow slopes I, II, and III of Fig. 6-3. 

In order to design the attenuator of Fig. 6-5, we need the 
equations of slopes I, II, and III, which are obtained as follows from 
Fig. 6-3. 
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FIGURE 6-5 

LABORATORY NON-UNEAR ATTENUATOR 
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Slope I : Eo = ^amp) ~ 29 

Slope II : Eq = ,354 E^^^^ - 9,91 

Slope III : Eo = .124 E(gj„p) ♦ 14.7 



Since the maximxam slope to be obtained from the attenuator shown in 
Fig. 6-5 is .5, we shall double the original value for the gain of 
elope I to compensate for this. Thus, the actual equations to be 
fitted hj the attenuator are 

Slope I : Eq - *5 £( aynp ) — 14*5 

Slope II : Eq = .177 E(amp) - 4.96 

Slope III ; Eo = *062 E(amp) 7.35 

The diode switching circuit shown in Fig. 6-5 will follow 
the double line curve shown in Fig. 6-7* This characteristic plus the 
negative 14.5 volts is then summed by the feedback amplifier to give 
the desired characteristic. In the actual circuit, all polarities are 
reversed since the DC amplifier involves a 180 degree phase change and 
only positive outputs of the DC amplifier are possible because of the 
unidirectional flow of current in the current source vacuvun tube. 

Specifically, the operation of the diode switching circuit is 
explained as follows. As the input, increases from zero along 

slope I, neither diode is conducting and the voltage at A fed into the 
summing amplifier is 



Ea 



% 

Rl + R4 



E(ainp) 



•5 ®(afflip) 



However, when A = .0175 rad as seen from Fig. 6-7, it is 
necessary to change this characteristic to follow slope II. Thus, at 
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o(votts) 




FIGURE 6-6 

DESIRED CHARACTERISTIC TO BE 
FOLLOWED BY NON-LINEAR ATTENUATOR 
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( volts) 




FIGURE 6-7 

ACTUAL ATTENUATOR CHARACTERISTICS 
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this point, V-1 goes into operation and the slope is now determined to 
a very close approximation by 



(R2 R3)(R4) 

Rl + (R2 * R3) (R/,^) 
R2 ♦ R3 + R4 



In the above expression, we have neglected the small resistance of the 
6 AL 5 diode (varying from 200 to 1000 ohms) and that of the biasing 
network. 



Going through elementary calculations it is found that when 

A = .0175 rad, the desired switching point, then the voltage at point A 

is - 7*3V. Thus, with the plate of V-1 set at - 7#3V as shown, the 

correct breakpoint is determined. Slope II is followed until A = .0335 rad. 

for the full ten degree operation. At this point, we desire V-2 to 

start conducting thus effectively shorting Rq. The slope for 

^ E(amp) 



is then quite closely found to be 



^2 % 

R2 % 

Rl + R2 R4 
R2 + R4 



.062 



Again, elementary calculations have shown that when A = .0335, 
point B of Fig. 6-5 is at - 9V. Thus, setting the plate at - 9V will 
give the proper switching action. 

Fig. 6-8 shows the response of the orthogonal vector summing 
system to 3 volt step functions without a non-linear attenuator in the 
feedback loop. A Sanborn pen recorder was used to plot the output of 
the system as a function of time in seconds. The gain was set so that 
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70o7 volts output gave the full ten degree operation. As explained 
previously, plot C of Fig. 6-4 shows the straight gain characteristic 
for full ten degree operation* Only at = 34 volts is the slope of 
the linear characteristic the same as that of the ideal curve. This 
is the point on Fig. 6-4 where d^ = d2 since the ideal curve is half 
way between A and B. 

From Fig. 6-8, it appears that the (DR)(ovss) referring to 
the linear second order system at Eq = 34V is about .6. 

To check this, we make use of the theoretical expression 



was calculated by putting a known small current through a torque 
generator winding and measuring the time taken to travel through a 
given angle. 



Thus, a rough check on the time it takes for the output to first pass 
through its final asymptotic value can be found by the expression. 




6-37 



Cd has been found to be 2.5 x 10^ at 200^F. This 

raci/ sec 



Substituting numerical values into 6-37 gives 

= .625 which is a good check on the experimental 




result. 



Wn has already been found to be 60 rad/sec for the ideal case. 



^p = i ^ ^ 



or 



7T 




6*“38 



In the case of the ideal system, Tp = ,033 sec, a value that 



98 



should correspond to the 34V point. Interpolating between the 30 V 
and 40V points on Fig. 6-8, this value is fairly well checked. 



Fig. 6-9 shows the same straight gain type of system as that 

of Fig. 6-8 except that Eo(max) corresponds to 5^ of shaft angle instead 

of the full 10^ in the previous case. Thus, at corresponding input 

levels, the damping ratio associated with the linear squared system of 

Fig. 6-9 is times those of Fig. 6-8. On the other hand, the un- 

\ 2 

damped natural frequencies of the levels of Fig, 6-9 are ^^2 times 
those of the corresponding levels of Fig. 6-8. A good checkpoint for 
comparing the two systems is at the 30V point. From Fig. 6-8, the 
damping ratio can be seen to be about ,7. On Fig. 6-9, it appears that 



the damping ratio of the same level is between and .5. But 



yfj 



times .7 is .49. Because of the small scale required of the recorder, 
the data is necessarily rough and thus it is quite proper to speak of 
any constants derived from the experimental results only in approximate 
terms. 



We shall now use the non-linear attenuator which was described 
earlier in this chapter. Figs. 6-lOA and 6-lOB show the responses of 
the system to a 3V step function for various operating levels. It is 
to be noted that for all levels from 3V to 50V, the dynamics are 
essentially constant. Even at the zero volt level, the response is 
greatly improved due to the very high gain at the very low levels. 

Maxima were limited to 50V in the experimental results since 
only one input was used and each input channel was designed to take a 
maximum of 50V. 

In conclusion, it is hoped that a fairly logical approach has 
been made in analyzing this particular non-linear system. Admittedly 
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FIGURE 6— lOB 
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the title, ”A Study of the Dynamics Associated With A Non-linear 
Closed Loop System,^’ left a certain amount of leeway in Just what 
was to be accomplished. However, the instrumentation of the ideal 
system in Chapter II in addition to the development of a method for 
designing a practical system shown in Chapter III should cover most 
of the background necessary to develop a practical system similar to 
that shown in Chapter VI. 

The author developed a firm appreciation for the difficulties 
involved in analyzing non-linear systems, and why they are avoided if 
possible. Although the use of perturbation theory is the most general 
method for attacking a non-linear system, each non-linear system is a 
separate problem in itself and must be treated as such. Thus, the 
method used to attack the Orthogonal Vector Summing System using 
microsyns probably will apply to very few if any, other non-linear 
systems • 

As can be seen from the data sheets, the results from the 
laboratory test system seemed to agree very well with the general 
theory developed prior to Chapter VI. Of course, there are probably 
better and easier methods for obtaining the required non-linear attenu- 
ator than the diode switching circuits. A light bulb could very well 
have the desired non-linear characteristics so that the pre-assigned 
tolerances on the dynamics will be maintained throughout the operating 
range, (except near zero). However, the diode switching circuit does 
not take much space and is fairly easy to design. Furthermore, by 
using more diode switches, the ideal square root curve can be approached 
as closely as desired. 
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APPENDIX 



Brief Explanation of Terminology 
Used in This Thesis 

Given: 1* An operating component, oco 

2* An input quantity, 

3* An output quantity, q^. 

Then the operating component performance function, (PP)(oc)» 
relates the input and output quantities as follows: 

% “ (Pf)(oc) 

can be separated into a static sensitivity, 

^(oc) ^qi»qoj » * constant of proportionality relating the input 

and output quantities, and a frequency function, (FP)(oc)* >'^ich tells 
how q^ is effected dynaaically by the operating component. 

For example, if we have a torque generator acting on a shaft, 
referred here €is a torque summing member, (tsm), than 

**(tsa) " ^(tg) [i^jm] ^^(in) • 

This equation is read, torque applied to the torque summing member is 
equal to the product of the sensitivity of the torque generator for input 
current squared and torque out multiplied by the square of the input 
current o 

The frequency function will be of the fora 

P“ ^ 

toto P“ b„.i ... 1 
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vrtiere p, as used in this thesis, is the Laplace transfora variable. 

For further explanation of the Laplace transform, the reader is referred 
to Reference k of the Bibliography, 
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